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Abstract. We study the hydro-dynamic limit for some conservative particle systems with 
degenerate rates, namely with nearest neighbor exchange rates which vanish for certain con- 
figurations. These models belong to the class of kinetically constrained lattice gases (KCLG) 
which have been introduced and intensively studied in physics literature as simple models 
for the liquid/glass transition. Due to the degeneracy of rates for KCLG there exists blocked 
configurations which do not evolve under the dynamics and in general the hyperplanes of 
configurations with a fixed number of particles can be decomposed into different irreducible 
sets. As a consequence, both the Entropy and Relative Entropy method cannot be straight- 
forwardly applied to prove the hydrodynamic limit. In particular, some care should be put 
when proving the One and Two block Lemmas which guarantee local convergence to equi- 
librium. We show that, for initial profiles smooth enough and bounded away from zero and 
one, the macroscopic density profile for our KCLG evolves under the diffusive time scaling 
according to the porous medium equation. Then we prove the same result for more general 
profiles for a slightly perturbed dynamics obtained by adding jumps of the Symmetric Sim- 
ple Exclusion. The role of the latter is to remove the degeneracy of rates and at the same 
time they are properly slowed down in order not to change the macroscopic behavior. The 
equilibrium fluctuations and the magnitude of the spectral gap for this perturbed model are 
also obtained. 



The purpose of this article is to define a conservative interacting particle system whose 
macroscopic density profile evolves according to the porous medium equation, namely the 
partial differential equation given by 



where A = Y2i<j<d®u an d m £ N\ {1}. This can be rewritten in the divergence form as 
d t p(t,u) = V(D(p~(t,u))V(p(t,u))) with diffusion coefficient D{p(t,u)) = m/o™" 1 ^,^). Note 
that D(p) goes to zero as p — > 0, thus the equation looses its parabolic character. 

One of the most important properties of the porous medium equation is that its solutions 
can be compactly supported at each fixed time or, in physical terms, that is has a finite speed 
of propagation. This is in strong contrast with the solutions of the classical heat equation. 
A non negative solution of the heat equation is always positive on its domain. A second 
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1. Introduction 



d t p{t,u) = Ap m (t,u) 
P(0,-) =Po(0 
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observation is that the solutions of the equation (jl.ip can be continuous on the domain of 
definition without being smooth at the boundary. The existence of this kind of solutions is a 
direct consequence of the degeneracy of D(p) as p — > 0. For a reference on the mathematical 
properties of the porous medium equation we refer to [11] and references therein. We also 
recall that such equation is relevant in different contexts in physical literature beyond the 
original motivation of describing the density of an ideal gas flowing isothermally through an 
homogeneous porous medium (corresponding to the choice m = 2). 

A microscopic derivation of the porous medium equation has been already obtained in [3] 
by considering a model in which the occupation number is a continuous variable. Here we 
study instead models with discrete occupation variables: our microscopic dynamics are given 
by stochastic lattice gases with hard core exclusion, namely systems of interacting particles on 
the d-dimensional discrete torus Tfj with the constraint that on each site there can be at most 
one particle. A configuration is therefore defined by giving for each site x E the occupation 
variable, rj{x) £ {0, 1}, which stands for empty or occupied sites, respectively. Evolution is 
then given by a continuous time Markov process during which the jump of a particle from a site 
x to a nearest neighbor site y occurs at rate c(x, y, rj). The choice c(x, y, rj) = 1 corresponds to 
the Symmetric Simple Exclusion process (SSEP) and, as is very well known, leads to the heat 
equation under diffusive re-scaling of time, namely D(p) = 1. In order to slow down the low 
density dynamics and obtain a diffusion coefficient which degenerates for p — > 0, we impose 
a local constraint (in addition to hard core exclusion) that must be satisfied in order for a 
particle jump to be allowed. This constraint is imposed at the level of rates: for any nearest 
neighbor couple (x, y) we fix the exchange rate c(x, y, rj) to be zero if rj does not satisfy a 
local constraint. The latter corresponds to requiring a minimal number of occupied sites in 
a proper neighborhood of (x,y). Since the typical number of particles in a given region is 
monotone with p, D(p) will decrease as p is decreased. At the same time the rates are chosen 
in order to satisfy the detailed balance condition with respect to Bernoulli product measure 
at any density (as for SSEP), namely the constraints do not add further interactions beyond 
hard core exclusion. The models we introduce belong to the class of kinetically constrained 
lattice gases (KCLG), which have been introduced and analyzed in physical literature since 
the late 1980's to model liquid/glass and more general jamming transitions (see for a review 
[9l HOj and references therein). In this context, the constraints are devised to mimic the fact 
that the motion of a particle in a dense medium (e.g. a molecule in a low temperature liquid) 
can be inhibited by the geometrical constraints induced by the surrounding particles 0. For 
most KCLG a degenerate diffusion coefficient is expected when p — ► but, for very restrictive 
choices of the constraints, the degeneracy could even occur at non trivial critical density. 

Here we provide the first derivation of the hydrodynamic limit for the simplest KCLG, 
the so called non cooperative KCLG. This means that the rates are such that it is possible 
to construct a proper finite group of particles, the mobile cluster, which has the following 
properties. There exists (at least one) allowed (i.e. with strictly positive rates) sequence 
of nearest neighbor jumps which allows to shift the mobile cluster to any other position. 
Furthermore, this allowed path should be deterministic, i.e. independent on the value of the 
occupation variables on the remaining sites. Finally, the jump of any other particle to a 
neighboring site should be allowed when the mobile cluster is brought in a proper position in 



Note that the role of particles and vacancy is usually exchanged in physical literature with respect to our 
convention: vacancies rather then particles are needed to facilitate motion. With this notation the diffusion 
coefficient degenerates at high rather than low density. 
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its vicinity. Therefore a configuration containing the mobile cluster can be connected to any 
other one with the same property by an allowed path. This, very loosely speaking, means that 
non-cooperative KCLG should behave like a re-scaled SSEP with the mobile clusters playing 
the role of single particles for SSEP. Thus their diffusion coefficient should degenerate only 
for p — > 0, when the density of mobile clusters goes to zero as power law with p. Indeed, in 
our proof of the latter result, the non cooperative property will play a key role since we will 
use it to provide paths which allow to perform particle exchanges. The use of similar path 
arguments for KCLG had already been exploited in [I], were the scaling with the lattice size of 
the spectral gap and log Sobolev constant for non cooperative models in contact with particle 
reservoirs at the boundary were derived. A similar case in which the diffusion coefficient does 
not vanish has already been studied in [8]. Finally, we stress that all along our proofs of 
hydrodynamics (both for the Entropy and Relative Entropy method) we use an additional 
property of the rates: the fact that they are of gradient type. A natural development of our 
work would be to generalize the present results to non cooperative KCLG with non gradient 
rates. 

An outline of the paper follows. In section[2]we introduce some notation, define our models 
and state the main results. The models we consider are either truly KCLG or perturbed 
models in which we add proper jumps of SSEP to remove the degeneracy of the exchange 
rates. In section [3] the proof of the hydrodynamic limit via the Relative Entropy method for 
the KCLG is presented. The proof of the hydrodynamic limit for the perturbed models via 
the Entropy method is described in section [H The proof of a Replacement Lemma needed 
for the Entropy method is postponed to the section [5j Finally, in section we study the 
spectral gap for both the unperturbed and perturbed dynamics. 



The models we consider are continuous time Markov processes r\t with space state Xd = 
{0, 1} T ^, where = {0,1,.., JV — l} d is the discrete d-dimensional torus. Let r\ denote a 
configuration in Xd , x a site in and rj(x) = 1 if there is a particle at site x, otherwise rj(x) = 
0. The elementary moves which occur during evolution correspond to jump of particles among 
nearest neighbors, x and y, occurring at a rate c(x, y, rj) which depends both on the couple 
(x, y) and on the value of the configuration rj in a finite neighborhood of x and y. Furthermore 
these rates are symmetric with respect of an x-y exchange c(x, y, rj) = c(y, x, rj) and are 
translation invariant. More precisely the dynamics is defined by means of an infinitesimal 
generator acting on local functions / : Xd ~ > ^ as 



2. Statement of results 




c(x,y,rj)r](x)(l - ri{y))(f(r] x > y ) - f{rj)) 



(2.1) 



\x-y\=l 



where \x - y\ = J2i<i<d \ x i 



yi\ is the sum norm in W 1 and 




(2.2) 



In the sequel we consider the rates 

c(x, x + ej,7]) = r](x 



ej) + rj(x + 2ej) 



(2.3) 
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where {ej,j = 1, -,d} denotes the canonical basis of W 1 and we will prove all the theorems 
for this choice. This, as we will prove, leads in the hydrodynamic limit to the porous medium 
equation (jl.ip for m = 2. Also, we can provide for any other m a proper choice of the rates 
such that all proofs can be readily extended leading in the diffusive re-scaling to the porous 
medium equation with the correspondent to. For instance in the case to = 3, the jump rates 
to be considered are 

c(x, x + ej, rj) = r]{x — ej)r](x + 2ej) + ij(x — 2ej)r/(x — ej) + r/(x + 2ej)r](x + 3ej). 

Note that both the choices of the jump rates taken above have the property of defining a 
gradient system, namely one for which the instantaneous current between the sites and 
e f 

W , ej (v) = <0,e j ,f,){ n {0)0.-f,(e j ))-f,(e j ){l-f,{0))} 
can be rewritten as a function minus its translation. This property will be a key ingredient 
when deriving the hydrodynamic limit. Also, both the models are non cooperative in any 
dimension according to the definition given in introduction. Consider for example the rates 
(12. 3p in one dimension. A possible choice for the mobile cluster is given by two particles at 
distance at most two. Let us describe the deterministic sequence of allowed moves (i.e. with 
strictly positive exchange rate) which we should perform to shift of one step to the right the 
mobile cluster when rj{x) = rj(x + e±) = 1, i.e. to transform rj into rj' with r]'{x + e\) = 
tj'(x + 2ei) = 1, rj'(x) = 7](x + 2ei) and rf (z) = rj(z) for z g" (x, x + e±, x + 2e±). First we 
make the move rj — ► r j x +e 1 ,x+2e 1 ^ which is allowed since c(x + e±,x + 2e\,rj) > rj{x) = 1. 
Then we perform the move f] x + e i> x + 2e i _> ^x+e 1 ,x+2e 1 ^x,x+e 1 w j 1 j c ] 1 j s a i so allowed, since 
c(x,x + ei , 7y x + e i ,3,+2ei > r jX+ei,x+2ei^ x _|_ ie{) = f](x + e\) = 1. The case in which we have 
instead the particles at distance two, r){x) = rj(x + 2e%) = 1, can be treated analogously. The 
second property which characterizes non cooperative models can also be readily checked: if 
we are given any two neighboring sites, (y, y + ei), the exchange of their occupation variables 
can be performed if we bring the mobile group of two particles in (y — 2e%,y — e%) since 
c (y>2/ + e iiV) ^ v(u ~ e i)- It is then possible to verify that any two configurations rj and 
rj' with the same number of particles, X] r /( x ) = ^2v'( x ) an d both containing at least two 
particles at distance at most two can be connected one to another via a sequence of allowed 
jumps. 

Let v a be the Bernoulli product measure in Xd with a G (0, 1). Since c(x, y, rj) = c(y, x, rf) 
Vx, y G T^, the measures v a are reversible for this process Va, as for SSEP. By the degeneracy 
of the rates, other invariant measures arise naturally. For example in the one dimensional 
setting, any configuration rj such that the distance between the position of two consecutive 
occupied sites is bigger than two has all the exchange rates which vanish. Therefore it is a 
blocked configuration and a Dirac measure supported on it is an invariant measure for this 
process. 

Let Ejv^ denote the hyperplane of configurations with k particles, namely 

z N ,k = {v e x% ■ E ^ = k ^ ( 2 - 4 ) 

XET% 

which is invariant under the dynamics. We say that O is an irreducible component of ^N,k 
if for every rj, £ G O it is possible to go from rj to £ by a sequence of allowed jumps. For 
SSEP, the hyperplanes ^n,k are irreducible components for any choice of k and N. In the 
presence of constraints, a more complicated decomposition in general arises. For example in 
d = 1 with the rates (j2.3|) . the above observation on blocked configurations and on the non 
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cooperative character of the model, leads to the following irreducible decomposition for the 
hyper planes. If k > N/3, Sjv^ is irreducible. Instead, if k < N/3, Sjv,fc is reducible and 
decomposable into the irreducible component which contains all configurations with at least 
one couple of particles at distance at most two plus many irreducible sets each containing 
only a blocked configuration 

The irreducible decomposition in dimension d > 1 is more complicated. In this case 
the model is still non cooperative and a possible mobile cluster is given by a d-dimensional 
hypercube of particles of linear size 2. For any choice of the spatial dimension d, it is possible to 
identify a constant C(d) < oo such that the hyperplane £jv,A; is irreducible for k > C(d) (N/3) d , 
while it is reducible in several components for k < C(d)(N/3) d . In this case we have: (i) the 
irreducible component which contain configurations with at least one d-dimensional hypercube 
of particles of linear size two plus all configurations that can be connected to these; (ii) 
irreducible components which contain single (blocked) configuration; (iii) other irreducible 
components which contain neither blocked configurations nor any configuration belonging to 
(i) . An example of irreducible set of the third kind for the rates (12. 3ft in d = 2 is for example 
the one that contains all configurations which have two particles at distance smaller or equal 
to two on a given line, for x = (xi,X2) such that rj{x + ei) = r](x) = 1 or rj(x + 2ei) = r](x) = 1 
and are completely empty V y = (2/1,2/2) which do not belong to the same line, namely yi / xi- 

In order to investigate the hydrodynamic limit, we need to settle some notation. Define 
the empirical measure by: 



Trf (du) = 7r N (r] t ,du) = — ^ rj tN a(x)6^(du), (2.5) 



where 5 U denotes the Dirac measure at u. 

Let T d denote the d-dimensional torus. Fix now, a initial profile po : T d — > [0, 1] and denote 
by (p n )n a sequence of probability measures on \d ■ 

Definition 1. A sequence (p n )n is associated to po, if for every continuous function 
H : T d — > M. and for every 6 > 



lim n N 



|^ E " j Td H(u)p ( U )du\ >d]=0. (2.6) 



Our goal consists in showing that, if at time t = the empirical measures are associated 
to some initial profile po, at the macroscopic time t they are associated to a profile pt which 
is the solution of the hydrodynamic equation (jl.ip . 

Fix e > and let p : T d -> [0, 1] be a profile of class C 2+e {T d ). By Theorem A2.4.1 of [5], 
equation (jl.ip admits a solution that we denote by p{t, •) which is of class C 1+e,2+e (IR+ x T d ). 

Here we also have to impose a bound condition on the initial profile, as the existence of a 
strictly positive constant 5$ such that 



Vn G T d , S < p (u) < 1 - 5 . (2.7) 

po(-) 

N 



Let 

v^iLtA be the product measure in \d suc h that: 



v 1 ; o (-){v,v(x) = l}=Po(x/N). 
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For two measures p and v in Xd denote by H{p/v) the relative entropy of p with respect 
to v, defined by: 

H{ji/v) =sup{ J f dp -log J e f dv}, (2.8) 
where the supreme is carried over all continuous functions. 

Theorem 2.1. Let p : T d — > [0,1] be a initial profile of class C 2+e (T d ) that satisfies the 
bound condition (2.1). Let (p N )^ be a sequence of probability measures on \d such that: 

H (» N /»M.)) = °( Nd )- ( 2 - 9 ) 

Then, for each t > 

7rf (du) > p(t,u)du 

N—>+oo 

in probability, where p(t, u) is a smooth solution of equation ( fi.il) . 

The porous medium equation (jl.ip presents an interesting behavior for profiles which vanish 
in some region of their domain. The previous method does not give us information in this 
situation since in order to apply it the initial profile has to satisfy condition (|2.7p . When these 
conditions are not verified, another possibility to derive the hydrodynamic limit is to use the 
Entropy method which is due to Guo, Papanicolau and Varadhan [3]. However the latter 
strongly relies on the ergodicity of the Markov process restricted to an hyperplane which, as 
we discussed above, does not hold for the process generated by Cp. One way to overcome this 
problem is to perturb slightly the dynamics in such a way that the frozen states are destroyed, 
ergodicity is restored and the macroscopic hydrodynamic behavior still evolves according to 

More precisely we consider a Markov process with generator given by 

Cf = C P + N e - 2 C s (2.10) 

where < 6 < 2, Cp was introduced in (12. ip and Cs is the generator of the SSEP, which acts 
on local functions / : Xd as 

(£sf)(ji)= E ^v(x)(i-v(y))(f(v x ' y )-f(ri)), 

\x-y\=l 

where rf' y was defined in (|2.2p . Since the Bernoulli product measures u a are invariant for both 
the processes generated by Cp and by Cs, they are also invariant for the perturbed process 
generated by Cg . Thanks to the restriction < 6 < 2, the generator of SSEP is slowed 
down and the Markov process generated by Cp and C^ will have the same hydrodynamic 
limit. Furthermore, since on any configuration each nearest neighbour exchange rate is strictly 
positive, the Markov process generated by C^ is ergodic on each hyperplane with fixed particle 
number and we can apply the Entropy method to derive its hydrodynamic limit. 

For a probability measure p on Xd ■> denote by W 6 ^ N = the probability measure on the 
space D([0,T],Xd ), induced by the Markov process with generator C^ speeded up by ./V 2 
and with initial measure p; and by the expectation with respect to P^. 

We start by introducing the definition of weak solutions of equation (jl.ip . 

Definition 2. Fix a bounded profile po : T d —> R. A bounded function p : [0, T] x T d — > M. 
is a weak solution of equation U.l\) if for every function H : [0, T] x T d -» R of class 
C 1 > 2 ([0,T]xT d ) 
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dt / du{p(t,u)d t H(t,u) + (p(t,u)) 2 dlM(t,u)} 
Jjd i<i<d 

+ / p (u)H(0,u)du = / p{T,u)H(T,u)du. (2.11) 

The Entropy method requires the uniqueness of a weak solution of the hydrodynamic 
equation. This is a consequence of Theorem A2.4.4 of [5] together with the fact that there is 
no more than a particle per site. 

Theorem 2.2. Let p$ : T d —* [0, 1] and (p N )n be a sequence of probability measures on Xd 
associated to the profile po- Then, for every t > 0, for every continuous function H : T d — * R 
and for every 5 > 0, 

1 E H^) m (x)- j^H(u)p(t jU )du 



lim ] 



> 5 



0, 



where p(t,u) is the unique weak solution of equation \l. 

Once we have established the Law of Large Numbers for the empirical measure for the 
process with generator Cff , the next step is to obtain the Central Limit Theorem starting 
from the invariant measure v p . 

For each z > (resp. z < 0) define h z : T d — > M by /i z (n) = \/2cos(27rz • u) (resp. 
h z {u) = \/2sin(27rz • u) and let /to = 1- Here • denotes the inner product of M. d . It is well 
known that the set {h z , z G is an orthonormal basis of L 2 (T d ). In this space consider the 
operator Q = 1 — A. A simple computation shows that Vlh z = "y z h z where "y z = 1 + 4vr 2 1 1 ^ 1 1 2 . 

For a positive integer k, denote by Tik the space obtained as the completion of C°°(T d ) 
endowed with the inner product defined by < /, g >k=< f, Q k 9 >■ Let 7i-k denote the dual 
of Tik with respect to the inner product of L 2 (T d ). 

Recall the definition of the empirical measure in (|2.5p where r\t denotes the Markov process 
with generator C e N , with < 9 < 2. As we want to investigate the fluctuations of this measure, 
fix p > and denote by 3^ the density fluctuation field acting on smooth functions H as 

,Af/m _ 1 —( x 



y?(H) = ^/2 E H(-)( Vt Mx)- P )- (2-12) 

Fix a positive integer k and denote by D([0,T],Ti.-k) (resp. C([0, T], Ti—k)) the space of 
Tl-k functions, that are right continuous with left limits (resp. continuous), endowed with the 
uniform weak topology. Denote by Qn the probability measure on D([0,T],7i-k) induced by 
y N and v p . 

Theorem 2.3. Fix an integer k > 3. Denote by Q be the probability measure on C([0, T], Tl-k) 
corresponding to the stationary generalized Ornstein- Uhlenbeck process with mean and co- 
variance given by 

E 3 {Y,(H)Y,(G)} = V ^">;' m I */ *5 W G W exp{-i^} 
y/8irp(t-s) Jro JRd 1 8{t-s)p) 

for every < s < t and H , G in Tik- Here H and G are periodic functions equal to H, G in 
Then, (Qn)n converges weakly to the probability measure Q. 
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The proof of this theorem is very close to the one presented for the Zero-Range process in 
[5 J and for this reason we have omitted it. We note that, since the proof of the Boltzmann- 
Gibbs in [5] relies on the ergodicity of the Markov process, we can use the ergodic properties 
of the generator Cs to obtain the result. We also remark if one considers the process rjt with 
generator Cp, we obtain in the limit the same Ornstein-Uhlenbeck process as described above. 
We only stress that, in this case, we derive the Boltzmann-Gibbs Priciple in a similar way as 
we do in the proof of the One Block estimate via the Relative Entropy method. 

Finally, in section [6] we investigate the spectral gap for the process in finite volume for 
d = 1. Our aim is to study the dependence on the size of the system of the spectral gap of 
the process. This, as is very well known, scales as 1/N 2 for SSEP on all hyperplanes £jv"fcj 
uniformly in k. For our models, due to presence of constraints, the uniformity in k is certainly 
lost, see remark \6. 2 1 

In order to illustrate our results we need to introduce a few additional notation. Fix an integer 
N and denote by A at the box of size N, A^r = {1,2... A^} and by \ N the configuration space 
X N = {0,1}^. In order to define the generator on \ N we could use definition 12.11 with the 
sum restricted to x, y G A^r. However, some care should be put when defining the jump rate 
for sites close to the boundary of A^v, since c(x,y,rj) as defined in (|2.3j) depend not only 
on the configuration on x and y but also on their neighbouring sites which can be outside 
Atv- We denote by dAjy the boundary set including all sites which do not belong to A^r and 
are nearest neighbour to at least one site in A^r, dAjy = {x G Z : x G" A^r, d(x, A^r) = 1}, 
where as usual the distance between a point x and the set A at is the infimum of the distances 
between y G A^r and x. A possible way to define the finite volume generator is to imagine 
that the configuration in the boundary set is frozen to a reference configuration, a, and to 
define the finite volume rates as c a (x,y,n) = c(x,y,r] ■ a) where c are the rates in (j2.3[) and 
rj ■ a G {0, i}I a jvI+I sa jv i s configuration equal to rj(x) on sites x G An and to cr(x) on 
x G <9Atv- In the following we make the choice a(x) = in x G dAjy and we denote by C>p,\ N 
and ^CgAjv ^ ne Markov processes with this choice corresponding to (|2.ip and (|2.10p . For sake 
of clarity, we explicitly write the second one 

(£%A N fKv) = E (c(a;,x+l,»7)+JV fl - 2 )77( a :)(l-T / (x+l))(/(Tf i ' a! + 1 )-/( J 7)) (2.13) 
xeAjv\{i,JV-i} 

+ £ (c(x,x - 1,7?) + A^- 2 )rKx)(l - V (x - IW^ 1 ) ~ f(v)) 
xeA N \{2,N} 

+ (r,(S) + N e - 2 ) (r/(l) - r/(2)) (/(r/ 1 ' 2 ) - f( V j) 

+ (r,(N - 2) + N 9 - 2 ) (rKiV - 1) - V (N)) (/(r/^ 1 ^) - /(r?)) 

where c(x,y,rj) and rj x ' y was defined in (12. 3h and (12. 2p . respectively. 

Let, with a slight abuse of notation, Sjv,fe denote again the hyperplanes with k particles, 
namely those in (12. 4h but with the sum running over Ajv- For each k, the Markov process 
generated by Cq \ n is irreducible on £jv,fc- The same holds for the process generated by £p,A N 
but only for k > N/3. This can be again proved by using the fact that the model is non 
cooperative with two particles at distance at most two being a mobile cluster. In both cases 
the unique invariant measure is the uniform measure, VN,k- 
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For a generator C with invariant measure fx, denote by \n(£) its spectral gap, defined by 

\ N {C) = mf — r, 
/gl 2 (m) Var(/,/x) 

where T>c(f,p) denotes the Dirichlet form defined by 

©r(/,M)=/ -f(v)£f(v)ridv)- (2-14) 

In the following we will also use the shortened notation Dp(f,p) and T>e(f,p) to denote the 
Dirichlet form with generator Cp,A N and £$\ N , respectively. Let p = k/N. We obtain that: 

Proposition 2.4. Fix k > AT/3. For the Markov process with generator Cp,^ N , there exists 
a constant C that does not depend on N nor k such that 

x n(£$,a n ) > Ajv(£p,Aiv) > C — !^ ■ 



Proposition 2.5. Fix k < N/3. For the Markov process with generator Cg^ N , where 6 = 1, 
there exists a constant C that does not depend on N nor k such that: 

3. The Relative Entropy Method 

In this section, we prove Theorem [231 Let p : T d -» [0, 1] be a profile of class C 2+e {T d ) 
that satisfies the bound condition f|2.Tj) and let {p n )n be a sequence of probability measures 
on Xd that satisfies (|2.9p . Fix a time t > 0. Denote by S^ ,p the semigroup associated to the 
generator Cp speeded up by iV 2 and by p^ the distribution of the process at time t. Denote 

Pit,-) 



by u N ( . x the product measure with slowly varying parameter associated to the profile p(t, •): 



^ it> . ) {v,v(x) = l}=p(t,x/N). 



It is well known that in order to prove the Theorem 12.11 it is enough to show: 

Theorem 3.1. Let po : T d — > K be an initial profile of class C 2+€ (T d ) that satisfies the 
bound condition \2. 7| ) and (p n )n a sequence of probability measures in \d that satisfies the 
condition Then, for every t > 

H(p?/^ (t , ) )=o(N d ), 

where p(t, u) is a smooth solution of 

Proof. Fix q £ (0, 1) and an invariant measure v a . Let 



dV N rl n 

= ex p{ E r,(x)\(t,x/N)} 



Since the measures v^ t ^ and v a are product, it is very simple to obtain an expression for 
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where 

f p(t,u)(l - a) 



^ , x , /P\hU)[l — a) \ 



a(l - p(t,u))J' 

and Zf is a renormalizing constant. 

In order to prove the result, we are going to show that 

ft 



H N (t) < o(N d ) + - I H N (s)ds 
7 Jo 



for some 7 > 0, and apply Gronwall inequality to conclude. 

There is a celebrated estimate for the entropy production due to Yau [12J: 

d t H N {t) < f _ dt log^( V )} f N (t)( V )v a (d V ), (3.1) 

where C* p is the adjoint operator of Cp in L 2 (u a ). 

Here and after, for a local function / we use the notation f(p) = E u [f(i])]. By simple 
computations together with the One-block estimate whose proof is presented at the end of 
this section, we can rewrite (ip^ (rj)) -1 {N 2 Cpip^ (77) — dtipj* (r])} as 

d 

Y.Y, 3 ^ X I N ) { W " *0><t> - h'{p(t, x/N))[ v l (x) - p(t, x/N)] } 

x£T d N i=l 



+ 



a 

Y, zZidu.X^x/N^^^x)) - ~g(p(t,x/N)) - ~g\p(t,x/N))[ V \x) - p(t,x/N)}}, 



x€T% 3=1 

plus a term of order o(N d ), where 

hj(v) = rifflMej) + viOH-ej) - r ? (-e J )7 7 (e j ), (3.2) 

g J ( V ) = c(0,e 3 , V )( V (0)-r ] (e j )) 2 (3.3) 

and 

Repeating standard arguments of the relative entropy method, the result follows. We refer 

the reader to chapter 6 of [5] for details. □ 

3.1. One-Block estimate. 



The main difficulty in the derivation of the One-Block estimate is the fact that, as already 
discussed, the process is not ergodic on hyperplanes with a fixed particle number. In order 
to overcome this problem, we separate the set of configurations into two sets: the irreducible 
component that contains all configuration with at least one d-dimensional hypercube of par- 
ticles of linear size 2 (and all the configurations that can be connected to them) and the 
remaining configurations. In the first case the standard proof is easily adapted, while for 
the second case we will use as a key ingredient the fact that this set has small measure with 



respect to u 



p(.t,-Y 
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Now we introduce some notation. Denote the also called Dirichlet form associated to the 
generator Cp and a measure fi in \d ■> by Dp{f,y) which is defined on positive functions by 

D P {f,fj)=J)p(y/J,f Jl ) (3.5) 

and S)p(/, fj) was defined in (|2.14|) . Let ff' P denote the Radon-Nikodym density of fj, N ' P (t) 

I r* 

t JO 



7 Jo ^ N S^ ,P ds with respect to u a . 



Lemma 3.2. (One-block Estimate) 

For every local function ip and for small 7 



limsuplimsup / V T x Vi^(n)f^ ,F \n)v a (dn) < [ H N (s)ds 



1 ' \y\<i 



where 



and with 77 (0) as defined in {3.J$ - 

Proof. Fix x G Tfj and denote by Q x \ the set of configurations in the box of center x and 
radius I containing at least one d-dimensional hypercube of linear size 2 which is completely 
filled: 

Q*,l = {v- E II 
yeC* z&Qy 

where Q y = {z : \zi - yi\ G {0, 1} Vi G {1, . . . , d}} and = {y : |y - x| < Z, Q y C T^}. 
We denote by S x i the irreducible set which contains Q x 1 (and all configurations that can be 
connected via an allowed path to one in £ x ,i) and we split the integral that appears in the 
statement of the Lemma into 

j W E ^ViMHe^m^ivMdr]) (3.7) 

xeT d N 

+ /iE T x V^{n)l {£ c } {n)f^ P {r,)v a {dn). (3.8) 



xdT% 

Note that the probability of the ergodic set £ x> i converges rapidly to one with /, indeed the 
following holds: 

^,.)(^)>i-(i-^y d (3.9) 

where we used the fact that the initial profile is bounded away from zero, namely p(s, it) > 60 
\/u G T d . Since H(/i N /v a ) = 0(N d ) and the entropy decreases with time, it holds that 
H(f t ' ) = 0(N d ). This implies that the Dirichlet form of f t ' is bounded above by CN d ~ 2 . 
For the term (|3.7p . since we are restricting on an irreducible set, we can repeat the standard 
arguments of the One-block estimate and conclude: 

limsuplimsup sup / — ^ V t x Vi^{j])1 {£x l} (r))f(rj)v a (dr)) = 0. 
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We now deal with the term (|3.8p and, to keep notation simple, we drop the integral with 
respect to time. The entropy inequality allows to bound it above by 



-)N d 



xer d N 



for every 7 > 0. 

The term on the left hand side is H^{s) /^N d . On the other hand, since is bounded 
and by Holder inequality, the second term of last expression can be bounded by 

N ~ ( 2l + IV* E lo S / exp {7(2/ + 1)^1^ ^,)}^.)^) 

x£T d N 

= nJi + iy £ log + im - 1) + 1). 

By using the upper bound on i/V ^ z ) which follows from (|3.9p and the inequality log(x + 
1) < x which holds true Vx, we can bound from above last expression by 

^ IF (exp{ 7 (2/ + l) d C] - 1)(1 - , 
which vanishes as I — > +00 provided 2^7 + log(l — S'q) < 0. □ 



4. The Entropy Method 

Now, we prove Theorem 12.21 The strategy of the proof is the same as given for the Zero- 
Range process in [5] . The main step is the derivation of the One-Block and the Two-Blocks 
estimate, which are presented in the fifth section. 

Fix a time T > 0. Let Jvi + be the space of finite positive measures on T d endowed with 
the weak topology. Consider a sequence of probability measures (Qn)n on D([0,T\,Ai+) 
corresponding to the Markov process ir^ as defined in (|2.5p . starting from /i . 

First we prove that (Qn)n is a tight sequence. Then we prove the uniqueness of a limit 
point, by showing that the limit points of (Qn)n are concentrated on trajectories of measures 
absolutely continuous with respect to the Lebesgue measure, equal to po(u)du at the initial 
time and whose density is concentrated on weak solutions of the hydrodynamic equation 
([Lip . By the uniqueness of these solutions we conclude that 7r^ has a unique limit point, 
concentrated on the trajectory with density p(t,u), where p(t,u) is the weak solution of 
equation (jl.ip . 

We divide the proof in several steps, to make the exposition clearer. Fix a smooth function 
H : T d —> K. Recall the definition of the empirical measure in (|2,5p and let 

By lemma Al.5.1 of |5j 

M^ H =< Tif , H > - < irg, H > - f N 2 C% < vrf , H > ds 

Jo 
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is a martingale with respect to the filtration Tt = o~(rj s ,s < t), whose quadratic variation is 
given by 

,t , d 



Z]^ EE W<*)*. (") 



where 



g ^(rj) = gj(rj) + (r/(0) - rj(ej)) 2 N e - 2 and ^(r?) as defined in (pT3|) . 

By elementary computations we can rewrite the integral part of the martingale as 

where 



di H( — )=N Z H( ^_^L ) + #( I Zl ) _ 2if, 

and ^(77) = + A^- 2 ? ? (0) and hj{rf) as defined in (fO) . 

4.1. Relative Compactness. 

Following the same arguments as for the Zero-range in [5] it is easy to show that (Qn)n is 
tight. 

4.2. Uniqueness of Limit Points. 

At first we note that all limit points Q of (Qn)n are concentrated on absolutely continuous 
measures with respect to the Lebesgue measure since there is at most one particle per site. 
The limit points are equal to po(u)du at the initial time by the hypothesis of (p n )n being 
associated to the profile p$. But, in order to prove that the limit points are concentrated 
on weak solutions of equation (jl.ip . we need to write the integral part of the martingale 
(|4.2p as a function of the empirical measure. This is the main difficulty in the proof of an 
hydrodynamical limit for a gradient system and we state it as a lemma: 

Lemma 4.1. (Replacement Lemma) 

For every 5 > and every local function tp 



lim sup lim sup ] 



/ J[d Yl T x V eNAVs)ds > 6 
J 1 



xeJ d N 



0, 



where Vi^ was defined in {37 

We postpone the proof to the following section in order to make the exposition more clear. 
To keep notation simple, in the previous statement and hereafter we write eN for [eN], it's 
integer part. 

Lemma 14 . 1 1 means that 

1 x - x - „/ X 



J2H(s,^){r x h j ( Vs )-h( V f(x))}ds 



vanishes in probability as N — > +00 and then as e — > 0, for every continuous function H. 
This is the replacement that permits to close the integral part of the martingale in terms of 
the empirical measure. For details we refer the reader to chapter 5 of [5]. 
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5. Replacement Lemma 

This section is devoted to the proof of Lemma 14.11 Before introducing the proof we give 
some fundamental inequalities that are used in the sequel. 

Fix a density a and an invariant measure u a , as reference. Let (/j, n )n be a sequence of 
probability measures on Xd and denote by the semigroup associated to the generator Cq . 
Denote by the density of /j, N with respect to v a . As the process is evolving on a torus 
we have that 

H(fi N /v a ) < C(a)N d , 

and since the entropy decreases with time we obtain that H(/j, N /v^) < C(a)N d , see 

Proposition Al.9.2 of [5]. Denote by (Cff)* the adjoint operator of in L 2 (u a ). Since the 

■A" 
t 

N _ »r2/ rN\* rN _ AT 2 pN f N 



process is reversible with respect to u a , / t satisfies the following equality: 



d t f»=N*(£$yf»=N^f« 



with initial condition — 



dug 



For / : Xd ~~ * reca U the definition of the Dirichlet form D(f,u a ) = D(^/J,u a ) and 
£>(/, Va) defined in (|2. 14[) . Since there exists a constant C, such that H^if^) < CN d , it can 
be proved that for all t > (see [5], section 5.2), that 

HN (lff" ds ) ^ CNd ' D e {^jj?ds,v a ) < (5.1) 

The proof of the Replacement Lemma relies on the well-known One-block and Two-blocks 
estimates. At first we reduce the dynamical problem, since the function depends on a trajec- 
tory, to a static one using the estimates obtained above. For that, we need to introduce some 
notation. 

Let fJL N (T) be the Cesaro mean of n N , namely: 

» N (T) = ±J% N S t N dt 

and fj? the Radon-Nikodym density of fi N (T) with respect to v a . Recall that in the beginning 
of this section we have obtained estimates for the entropy and the Dirichlet form of f^. As 
a consequence, to prove the Replacement Lemma it is enough to show that 

limsuplimsup sup (^^r x V eN ^{rf)f(ri)v a (dr])=Q, (5.2) 

£ ^0 N->+co Dg(f,u a )<CN d - 2 J ^ 

for every C < oo. This is a consequence of the following two results: 

Lemma 5.1. (One-block Estimate) 
For every finite constant C , 

limsuplimsup sup / i V T x Vi^(rj)f(r))v a (drj) = 0, 



where Vi^u was defined in \3. 6\) . 
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Lemma 5.2. (Two-blocks Estimate) 
For every finite constant C, 



limsup limsup limsup sup sup 

l^+oo e^O N^+oo D g {f,u a )<CN d - 2 \y\<eN 

[ W E Irffr + V)- V eN (x)\f(v)Mdri) = 0. (5.3) 

5.1. One-Block Estimate. 

Since in this case we are considering the perturbed process Cg , we can use the One- 
Block estimate for the SSEP speeded up by iV , which is enough to conclude the One-Block 
estimate for the perturbed process, the idea to proceed is the following. By the definition of 
the Dirichlet form (|2.14p and the generator : 

D e (f, u a ) = D P (f, u a ) + N e - 2 D s (f, v a ), (5.4) 

where Dp, (resp. D$) denotes the Dirichlet form as defined in (13. 5ft , Since the Dirichlet form 
is always positive, and we are restricted to densities / for which Dq(/, u a ) < CN d ~ 2 , it holds 
that: 

£>s(/,0 < N 2 - e D e {f,u a ) < CN d - e 
Following the same arguments as for the Zero-Range in section 5 of [5 J , it is not hard to show 
that the modified One-Block estimate 

limsup limsup sup / — ^ V] T xVi )1 p{i])f(ri)v a (dri)=0. 

,^ +00 7V^ +O o D s (f,u a )<CNd-B J ^ 

holds for the SSEP. The main difference between the proofs, comes from the bounds on the 
Dirichlet forms. Having the bound Ds(f, v a ) < CN d ~ 2 it provides the estimate D l s (fi, u a ) < 
C/N 2 , where // is the conditional expectation of / with respect to u-algebra generated by 
{77(3;), | x| < /}, while the bound Ds(f,u a ) < CN d ~ e , provides the estimate D l s (fi,u a ) < 
C/N 9 , which is enough to conclude the standard proof of the One-block estimate, since 
9 > 0. 

5.2. Two-Blocks Estimate. 

In order to show (|5.3p it is enough to prove that 

lim sup lim sup lim sup sup sup 

l^+oo iV^+oo D e (f,v a )<CN d - 2 2l<\y\<eN 

J ' Jfd E \v l (x)-r] l (x + y)\f(r])u a (dr])=0. 
We can write this integral as 

V(o)-^(y)l/faK(<fy), 



where / denotes the average of all space translations of /. 

In this case, we are able to separate the integral over configurations with small density of 
particles from the ones with high density. The first case is easily treated. For remaining, since 
we are taking configurations with high density of particles, we are able to construct a coupled 
process, in which each marginal evolves according to the SSEP and they are connected by 
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jumps of Cp. The main difference between the proof for the Zero-Range Process in [5] comes 
from the estimate in the Dirichlet form of this process, which in this case we are able to show 
that is bounded by C(d, l)e 6 and is enough to conclude. We divide the proof in several steps, 
in order to make the exposition clearer. 

5.2.1. Cut off of small densities. 

At first we note that if we restrict last integral to the small density set 

\x\<l \x— y\<l 

it is bounded above by ^2i+\) d ' wn i cn vanishes as I —>■ +oo. So, in fact we just have to consider 
the integral over the complementary set 11^. 

5.2.2. Reduction to microscopic cubes. 

Here we need to introduce some notation. Fix a positive integer I and a site x, denote by 
Ai(x) the box centered at x with radius /. Let x 2 ' 1 denote the configuration space {0, l} Ai (°) x 
{0, l} Ai (°), £ = (£i, £2) the configurations of x 2,1 and by u 2,1 the product measure u a restricted 
to x 2,1 - Denote by f y j the conditional expectation of / with respect to the cr-algebra generated 
by {r]{z),z G Ai(0,y)}, where 

A { (o,y) = Aj(0) u A { (y). 

Since, 77' (0) and rf(y) depend on t](x), for x G A/(0,y), we are able to replace / by f V) i and 
rewrite our desired limit as 

limsup limsup limsup sup sup 

Z^+oo e^O N-^+oo D e (f,is a )<CN d - 2 2K\y\<2eN 

\c[(o)-m\fyM> 2 a l m = o. 

5.2.3. Estimates on the Dirichlet form. 

At first we need to introduce some notation since we are working in the space J)q y . Recall 
from (|2.14p the definition of the Dirichlet form of /. Since the Bernoulli product measures 
are homogeneous, we can define the Dirichlet form of / as: 

D e {f,V a )= ^M'Va), 

x,zeJ% 
\x—z\=l 



where 



<*(/, = IxM, "«) + N e - 2 I s XiZ {f, i/ a ), (5.5) 
l xAf^o) = \ c(x,z,7]) \J f{rf' z ) -y/f(rj) v a (dn) 



and 



By the definition of / and since the Dirichlet form is convex, this implies that Dg(f, v a ) < 
Do(f,v a ). 

The main step in the proof, consists in obtaining an upper bound for the Dirichlet form of 
f y j from the upper bound of the Dirichlet form of /, in such a way that we can use the ergodic 
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properties of the Markov process. The idea is to obtain a limit density with Dirichlet form in 
A/(0,y) equal to 0, and then by using the irreducibility we can decompose that density along 
the hyperplanes. As the boxes Aj(0) and A[(y) have no communication, we cannot define the 
Dirichlet form in Aj(0, y) as the sum of the Dirichlet forms in Aj(0) and in A/(y), we must 
add a term that connects both boxes. Define on positive densities / : x 2,1 ~~ * 



(5.6) 



x,2GA((0) 
|cc-«|=l 



i,jzeAi(o) 

|X-2| = 1 



where 



2. l\ 

a J 



r-Hf v 

± x,z\J 1 U t 
j2,l , f u 2,l 



/(#*,&)- sfmM i&idZi,^), 



0,y 



/ 2 / 

We still have to define the term that connects both boxes, namely Iqq(/, v& ). Notice that 



n 



0,y 



n c un c v , where 



Q = {r ] e{0,l} Al{0) : £ ^)<2 d + l} 



and 

Then we define 



{r?e{0,l} Ai ^ : £ ^ < 2d + !}- 

zeA;(j/) 



j = l 



w i=i 



where £^ ,=l: = £j db <9o and <9o is the configuration with one particle at site and in the rest 
empty. 

This Dirichlet form corresponds to a particle system on Aj(0) x Aj(0), where the marginal 
processes evolve as SSEP and where particles can jump from the origin of one marginal process 
to the origin of the other and vice- versa, according to the jumps of the generator Cp. Using 
Schwarz inequality and the definition of f y i, we have that: 

Ix?z(fy,l, v 2 J) < I% z (/, v 2 J) and Ixfz(fy,l, "¥) < ^(/, v* 1 ) 

Then, by equality (|5.5|) and the estimate on the Dirichlet form of /, it holds that Dg{f , v a ) — 
CN d ~ e which together with l/N < e, implies that: 



x,zGAi(0) 
\x-z\=l 



i,zGAi(0) 
\x—z\=l 



< C(d,l)N~ 



(5.7) 
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In last expression C(d,l) is a constant that depends on d and In what follows it may 
vary from line to line. 

It remains to obtain an upper bound for lbo{f y ,i)- By using Schwarz inequality and the 

— / — 2 / 

definition of f y> i, we can bound Iq (fy t i, vd ) by 

d 



o j=l 



/(??) i/ a (d7/) 



+ 



„ d 



y + e^r?) V f(v°> y ) ~ V /fa) 



Let Eq-' (f,u a ) j = 1,2 denote the first (respectively second) expectation above. In order to 



keep the proof clear we are going to estimate E l ^(f,v a ) and state it as a lemma. We note 
that a similar argument provides the same bound for E l Q 2 y (f, u a ). 

Lemma 5.3. Let f be a density such that Dg(f,u a ) < CN d ~ 2 and let y £ : 21 < \y\ < 
2eN. Then, 



E^ y (f,v a )<C(d,l) 

Proof. The idea consists in expressing the exchange rf ,y by means of allowed nearest neighbor 
exchanges rj x,x+l of the generator Co . Since the integral is restricted to we have for certain 
2 d + 1 particles in Aj(0). We discuss for definiteness the path in the case ry(0) = 1, r](y) = 0, 
the other possibility can be treated analogously. In order to bring the particle from to y we 
first move 2 d of the particles in A;(0) close to by means of the jumps that corresponds to 
Cs- Then we arrange them in order to form a d-dimensional hypercube of linear size 2, which 
is a mobile cluster. Now we can shift this cluster plus the particle originally in in each of 
the d directions by using only the jumps in Cp and we bring them close to y. In figured] and 
[2] we show as an example the path which allows to shift in the e\ and e2 direction the mobile 
2x2 square of particles (black circles) plus the particle originally in (grey circle). 



e . 



Figure 1. Moving in the direction e\ 



• * « 



• • i. 



1 



* 

* • 



Figure 2. Moving in the direction e2 
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Then we drop particle that was originally in in site y and bring back the mobile cluster 
alone to A;(0) again using moves in Cp. Finally, using the jumps of Cs, we put them in their 
initial positions. 

Note that the above path is uniquely defined only after choosing the 2 d particles in the box 
Aj(0) which are used to form the mobile cluster and whose existence in guaranteed by the 
fact that the integral is restricted to S7q. This will bring an entropy factor which corresponds 
to the number of possible initial positions for the 2 d particles, which is bounded by l d2d . 

Once for a given r/ the initial position of the 2 d particles which will be used to form the 
mobile cluster is fixed, namely x±, ..,x 2 <i, we let r 7 o,j/r 7 o,s _i...T\(rj) be the sequence of nearest 
neighboring exchanges that represents the above described path. We can therefore rewrite, 
for any function g, the square difference (g(r] 0,y ) — g(r])) 2 as a telescopic sum: 



v o,a 



fc-1 



fc-1 



fc=l 



where 7 ' y denotes the number of steps of the path which is less than C(d)(l + eN). 

Now we separate the part of the path in which the moves are performed via jumps in Cs 
from the one that can be performed by using Cp and, by using the elementary inequality 
(x + y) 2 < 2x 2 + 2y 2 , we bound last expression by 



C(d)l fc-1 fc-1 C{d)Ne fc-1 fc-1 

[E 9{r k l[n( v )) -g(Hn( v ))] +[ £ g(r k \{ 7?fa)) - g{ J] n( V ) 

k=C(d)l+l 



fc=l 



8=1 



i=l 



i=l 



i=l 



Note that in the previous expression we had to take into account the number of steps in the 
paths, which is order I for the terms involving moves of Cs and of order Ne for those using 
Cp, thanks to our choice of the path. Also note that the term on the left hand side of last 
expression, refers to all the jumps inside the box of size I: when bringing the 2 d particles 
together and when putting them in their initial position. 

We now apply the above inequality with the choice g = \fj and use again the Cauchy- 
Schwarz inequality on each single term and the fact that the rates of Cp (|2.3p are bounded 
from below by 1 thanks to our choice of the path for all the terms in the second telescopic 
sum. This leads to the bound 



<C{d)l £ 

aieAi(O) 
a 2 de'X;(0) 



1 



{a 1 =x 1 , 



f(r]) v a {dri) 



-C(d)Ne £ 

aiGAi(O) 



L {«i 



i 2 d 



£c(ei,ry) yj f(rf 



where {ej}j denotes the bonds that we use inside A;(0) when taking the 2 d particles whose 
initial positions close to the particle at the site 0, by the jumps of the exclusion 

process, while corresponds to the bonds used by the generator Cp when performing the 
rest of the path and c(e«, 77) denotes the corresponding jump rate. Since there are l d2,d chances 
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for the initial positions Xi, . . . 



x. 



,d, we can bound last expression by: 



l d2d C(d)l Yl ll x+ xtf,v«) + l d2d C{d)eN 4 P ,+i(/,^)- 



x,x+l€Ai(0) x,x+leAi(0) 

By the equality (|5.5|) and the bound on the Dirichlet form of /, it holds that: 



which together with Z < eA^, ends the proof. 



□ 



2 / 2 Z 

By the definition of the Dirichlet form D g ' (f,vd ) in (|5.6p . together with (|5.7p and the 
previous Lemma, we can restrict to densities / that satisfy 



D 2 /(f,u 2 J)<C(d,l)e e . 



So, to conclude the proof of the Two-blocks estimate it is enough to show that 



hmsuplimsup sup f |£(0) - £ 2 (0)|/(£)^(^) = 0, 



where the supreme is carried over densities with respect to v a ' ■ Note that the parameter e 
is appearing only on the bound of the Dirichlet form. The proof now follows the same lines 
as in the case of the Zero- Range process in [5], by starting to decompose the Dirichlet form 
along the hyperplanes and applying the equivalence of ensembles. For the detailed strategy 
see section 5.5 of [5]. 



In this section we analyze the magnitude of the spectral gap for the one dimensional 
generators on finite boxes, Cp,A N and Cg^ N , which have been defined in section 2. Note that 
the results below on the scaling of the spectral gap with the lattice size have not been used 
in the previous sections to derive the hydrodynamic limit. This was possible thanks to the 
fact that we were considering a gradient choice of the rates. The following analysis of the 
spectral gap can be regarded as a first step towards the analysis of the non gradient version 
of our models, e.g. for the choice 



which is a non gradient version with the same kinetic constraints as our original choice (|2.3p . 
namely c'(x, x + e^, 77) = if and only if c(x, x + ej, 77) = 0. 

For simplicity in the following we drop A^v from our notation. The main ingredient of our 
proofs will be a comparison with the spectral gap \(£lr) for the long range exclusion process 
and the use of path arguments. Let us start by recalling the definition of Clr and the result 
in [7] for its spectral gap. The action of Clr on local functions / : Ejvfc — ► K is given by 



6. Spectral Gap 





(6.1) 



z.j/sAjv 



where rj x ' y as defined in (|2.2|) . 
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Consider the Dirichlet form T>lr associated to £l_r, with respect to the uniform measure 
un k, given explicitly by: 



VLR(f,m,k) = -^ £ / (f(V x ' V )-f(ri)) 2 VNAdv)- (6-2) 

x,yeA N Js N,h 

Quastel in [7], obtained that 



Var(/, u Ntk ) < T>lrU, VN,k), (6.3) 

by computing precisely the eigenvalues of the generator Clr- In order to prove our results 
(Proposition 12, 4| and I2.5[ ) for the spectral gap of the process with generator C^ Vfc and with 
generator Cp for k > 1/3, we proceed in two steps. 

Fix an integer k such that the density is restricted to p = > |. With this restriction, 
for each n S Sjv,fe, there exists a couple of particles whose distance is smaller or equal to two 
and £jv,& 1S irreducible. Furthermore we will show that there exists a constant C that does 
not depend on N nor k such that 

VlrU, y N ,k) < — ^riV 2 2)p(/, u Njk ), (6.4) 
P- 1/3 

where T>p(f, i*N,k) is the Dirichlet that corresponds to Cp. 

Last result together with (|6.3p allows to conclude the following Poincare inequality: 



Proposition 6.1. Fix k > N/3. For the Markov process with generator £p,N a nd for every 
f G L 2 (vn^), there exists a constant C not depending on N nor k such that: 

Var(f,u N>h ) < -^—N 2 T> P (f,u N , k ). 

The second inequality in the result of Proposition 12.41 is an immediate consequence of 
last result. The first inequality follows from the fact that for any function / it holds that 

The case in which p = — < | is more demanding. Since the hyperplanes £jv,fc are no longer 
irreducible for Cp, the corresponding spectral gap is zero. A natural issue is determining the 
magnitude of the spectral gap for Cp on the restricted irreducible set of configurations with 
k particles and at least one couple of particles at distance at most two. In this case the 
invariant measure is no longer the same one as for the long range jumps and the comparison 
with the latter process is no more useful. Instead, we will here consider for k/N < 1/3 the 
spectral gap for the modified generator Cff which, thanks to the addition of the exclusion 
part, is ergodic on the hyperplanes Sjy.fc for any k and reversible with respect to VN,k- By 
comparing the Dirichlet form of C^ with the one of Clr, we show that: 

Proposition 6.2. Fix k < N/3. For the Markov process with generator given by Co , with 
6 = 1 and for every f £ L 2 (i^/v,fc)> there exists a constant C not depending on N nor k such 
that: 

Var{f,u Ntk ) < ^N® s (f,VN,k) + -iV 2 2) P (/, u N>k ). 
P 2 P 



Proposition 12.51 is an immediate consequence of Proposition [67 
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6.1. Proof of Proposition [67TT , 

Fix an integer k and suppose that p = > 4. The Dirichlet form associated to Cp is 
given explicitly by 

\x-y\=l 

where c(x,y,rj) and rj x ' y were defined on (|2.3|) and (|2.2|) . respectively. We have seen above 
that it is enough to show (|6.4|) . 

The idea consists in expressing the exchange rj x ' y , for each x,y £ An by means of a sequence 
of allowed jumps of Cp. Consider for example the case x < y and rj(x) = 1, rj(y) = 0. The 
restriction p > \ guarantees that there exists a couple of sites (a, b) both in A at such that 
77(a) = 77(6) = 1 and \a — b\ < 2. Once the couple has been chosen we shift it close to sites 
x — 2, x — 1. We have thus reached a configuration with r/(x — 2) = rj(x — 1) = i](x) = 1, and we 
now shift the three particles to (y — 3, y — 2, y — 1). The allowed path which performs the shift 
of one step to the right is the composition of the following three basic steps: 77 — ^ 771 = >q x,xJr \ 
771 — > 772 = rj x ~ 1,x and 772 — > 173 = rj :r ~ 2,x ~ 1 ■ When after a proper number of one step shifts 
we reach the configuration r] n with r](y — 3) = r](y — 2) = 77(77 — 1) = 1, we can perform the 
exchange r\ n — ► 7? n +i = ,?/ which corresponds to dropping the particle at y. Then we make 
a similar backward path, this time shifting the two particles plus the vacancy untill bringing 
the vacancy in x and finally we bring back the two particles to their original position (a, b). 

As we did for Lemma 15.31 we have to take into account the entropic term corresponding 
to the possible initial positions of the couple of particles (a, b). Some care is required at this 
point, since using the rough bound N for this original position would lead to an additional 
factor N in the inequality 16.41 Let us start by presenting a lower bound for the number of 
couples of particles at distance at most two 

Lower bound for the number of couples for k > A^/3 

Fix a configuration r\ £ £jv,jfc- Let B be the set of sites which are occupied and such 
that for each of them there is a particle at distance at most two, B = {x : r](x) = 1,77(3; + 
1) + r/(x — 1) + 77(2; + 2) + 77(3; — 2) > 1}. Let also A denote the remaining occupied sites 
A = {x : 77(2) = l,x £ Ajy \ B}. The following holds: 

k = \A\ + \B\ 

N > 3(\A\ — 1) + 1 + \B\, 

where for a set 5, |5| denotes the cardinality of the set. 

The first equality is obvious. In order to establish the second property we use the fact 
that, if x and y belong to ^4, then \x — y\ > 3. The inequality follows by organizing the k 
particles in the closest configuration which does not change the value of \A\ and \B\. The 
minimum number of couples of particles at distance at most two and such that there are not 
common sites among different couples is \B\/2, see e.g. the figure below in which there are 8 
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neighboring particles but just 4 couples. This, together with the above established relations 
among k, N, \B\, gives the following lower bound for the number of couples 

r]{z)ri{z + 1) + r)(z)r)(z + 2) > \B\/2 > -(k - N/3 - 2/3). (6.5) 

zeAjv 

Introducing this term in the Dirichlet form of the long range exclusion (16. 2D we bound it 
by: 

jm E E EE^+wv/wM- ( 6 - 6 ) 

^.S/GAjv V& s N,k zeAjv 1=1 

Let us now consider a term in the sum above 

2 



V(z)v(z + l)[f(ri x ' y )-f(v) 

and let x — {z + 1) = n and y — x = m. By using the construction sketched above for a 
the possible path which connects 77 to 77 x,y via allowed elementary exchanges, it is possible to 
define a sequence r\i for 1 < i < 7^ ' y with 7^ ,2/ = 5(m — 1) + 4(n — 1) + 2 with the following 

properties: r/i = 77, r/ 7 j,» = 7/^ and Vi there exists ij 6 Ajy such that r\i = ^_W ,a! W +1 an d 
c(x(i), r?j-i) > 0, namely the exchanges are permitted for the generator Cp. Therefore 

we can rewrite each r](z)r](z + l){f{r] x ' v ) — fiv)) 2 as the- telescopic sum: 

7?' w -l 

r ? (z)r ? (z + /)(/(r / ^)-/(r ? )) 2 = r / (z)r / (z + /)( £ f(Vn)-f(Vn-l] 

i=l 

By using this equality together with c(x(i),x(i) + 1, 77i_i) > 0, the fact that the size of the 
path t£ is of O(N) and applying Cauchy-Schwarz inequality we can finally bound (|6.6p from 
above by: 

a n 2 

m\ E E E E tf*M* + E * * + M/fo* ) - 

x,j/GAjv V^N,k z£A N 1=1 e, 

where 2j = {i, i + 1} denotes one bond that we have used when performing the path that 
takes the particle from x to y using the couple at the sites z and z + 1 and C is a constant 
independent on JV and /c. Estimating over all x and y G A^r, last expression is bounded above 
by 

C AT 3 

E E ^)£c(*,*+ mx/^ 1 ) - mfvsM- 

By the estimate ()6.5[) and since the number of particles is fc, last expression is bounded above 
by 
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6.2. Proof of Proposition 16.21 




P P 



(6.7) 



which is enough to conclude. As before, the idea consists in expressing a path from x to y 
using the admissible jumps of C$ . 

For this choice of k we are no more guaranteed that there exist two particles at distance at 
most two, which was a key ingredient to construct the path in the high density regime. The 
idea will be to make use of the simple exclusion jumps to construct such mobile clusters and 
then proceed as before in a similar way as we did in in the proof of Lemma 15.31 

Fix a distance j and denote by Bj the set of particles at distance at most j, Bj = {x : 
r](x) = 1,31 £ (— j, — 1) U (l,j) s.t. rj(x + I) = 1} and by Aj the remaining particles. Then, 
the following holds: 



This inequality is obtained in the same manner as before considering that now the minimum 
distance that the „4j-particles have to be is j + 1. 

As before, the minimum number of couples that one can have is \Bj\/2. By simple compu- 
tations we obtain the lower bound 



Introducing this inequality in the Dirichlet form of the long range exclusion (|6.2p , we bound 
it from above by: 



For each configuration r/ and each choice x, y, z and I, we can now construct a path which 
first bring together the two particles in z and z + 1 by using the jumps of the simple exclusion 
and then uses this mobile cluster to perform the exchange of occupation variables in x and y, 
as in the proof of Lemma 15.31 Here Cs, also denotes the generator of the Symmetric Simple 
Exclusion process restricted to the box A at, given on local functions by 



Since the jumps of the exclusion are used just to put the neighboring particles at a distance 
equal to two, the size of the path for this process is 0(j). With this purpose, write f(r] x ' y ) — 
f{rj) as a telescopic sum, use the elementary inequality (x + y) 2 < 2x 2 + 2y 2 and then the 
Cauchy-Schwarz inequality to bound last expression by: 



k = \Aj\ + \Bj\ 
N> (j + l)(\Aj\ -\) + \ + \Bj 






(£sf)(v)= £ lr,(x)(l-r,(y))(f(r, x >y)-m). 



\x-y\=l 



Cj 



J 




(6.9) 
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(IN - 

+ ^\b-\ E E E E + o E < + 1, /(^) - /(*7))V*(^ (e-io) 

where ej denotes the bonds that we have used when bringing the neighboring particles at a 
distance equal to two, while e« = {i, i + 1} denotes the bonds we have used when performing 
the remaining part of the path that takes the particle from x to y. Note, that there is a factor 
j multiplying the first expression which comes from the size of the path for the jumps of the 
exclusion, while for the other process the size of the path is of 0(N). 

First we deal with the jumps that concerns Cp. As before, we bound (|6.10p from above by 

E E E x + X > *?X/(^ +1 ) - /fa))%r,*fa)> 



and since r\ S E/v,fc, we obtain the bound: 



P ~ 1/j 

Now, we bound (|6,9p . by 



E E E^)E(/(^ +1 )-/^)V fc w- 



MB,, 

J V^N,k z6Ajv Z=l x 

Since the number of particles is k, we can bound last expression by: 

C \ P ,., N® s (f,v N>k ). 
(p- l /J> 

Reorganizing these facts together we obtain that: 
Optimizing over j, (|6.7|) follows. 

Remark 6.1. For sake of simplicity we have presented the spectral gap results only in the one 
dimensional setting. By a proper modification of the path arguments and an accurate estimate 
on the minimal number of mobile clusters it is possible to obtain for d > 1 an analogous result 
as the one in Proposition \2.J\ if the density p = k/N is such that k > C(d)(N/3) d . 



Remark 6.2. Taking for instance the density fluctuations field as defined in \2.12}) and the 

reference measure the Bernoulli product measure v p , by simple computations we obtain that 

Var(y t N (H),v p ) = p(l-p)\\H\\l 
while the Dirichlet form corresponding to Cp equals to: 

Vp(y?{H) 1 v p ) = ±- 2 p\l-p)\\H>\\l 

So, if we consider H , such that \\H\\?, = WH'W?,, then: 

N 2 

Var(y t N (H),v p ) = —® P {y?(H),v p ), 
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which implies that the spectral gap )^n(£p,a n ) < T^r- 27ns is in agreement with the bound 
that we have obtained in ^6. 7\ ), when considering the spectral gap with respect to the uniform 
measure. 

References 

[I] Bertini, L.; Toninelli, C. (2004): Exclusion Processes with Degenerate Rates: Convergence to Equilibrium 
and Tagged Particle. Journal of Statistical Physics 117, 549-580. 

[2] Ekhaus, M.; T.Seppalainen, T. (1996): Stochastic dynamics macroscopically governed by the porous 
medium equation for isothermal flow Ann. Acad. Sci. Fenn., Math. 21 309-352. 

[3] Feng, S.; Iscoe, I.; Seppalainen (1997): A microscopic mechanims for the porous medium equation. 
Stochastic Processes and their Applications 66, 147-182. 

[4] Guo, M.Z.; Papanicolau, G.C.; Varadhan, S.R.S. (1988): Nonlinear diffusion limit for a system with 
nearest neighbor interactions. Commun. Math. Phys. 118, 31-59. 

[5] Kipnis,C, Landim,C.(1999): Scaling Limits of Interacting Particle Systems, Springer- Verlag, New York. 

[6] Liggett, T. (1985): Interacting Particle Systems, Springer- Verlag, New York. 

[7] Quastel, J. (1992): Diffusion of Color in the Simple Exclusion Process. Communications on Pure and 

Apllied Mathematics XLV, 623-679. 
[8] Spohn, H. (1991): Large scale Dynamics of Interacting Particles. Springer- Verlag. 

[9] Ritort, F.; Sollich, P.(2003): Glassy dynamics of kinetically constrained models, Adv. in Phys 52 219-342. 
[10] Toninelli, C; Biroli, G.: Jamming Percolation and Glassy Dynamics to appear in J. Stat. Phys. - published 
online, ISSN 1572-9613 

[II] Vazquez, J.L. (1992): An introduction to the mathematical theory of the porous medium equation, in: 
Delfour, M.C. and Sabidussi, G. editors, Shape Optimization and Free Boundaries, Kluwer, Dordrecht 
261-286. 

[12] Yau, Horng-Tzer (1991): Relative Entropy and Hydrodyamics of Ginzburg-Landau Models, Letters in 
Mathematical Physics, 22, 63-80. 

IMPA, Estrada Dona Castorina 110, CEP 22460 Rio de Janeiro, Brasil 
E-mail address: patg@impa.br 

IMPA, Estrada Dona Castorina 110, CEP 22460-320 Rio de Janeiro, Brasil, CNRS UMR 6085, 
Universite de Rouen, UMR 6085, Avenue de l'Universite, BP. 12, Technopole du Madrillet, 
F76801 Saint-Etienne-du-Rouvray, France 

E-mail address: landim@impa.br 

Laboratoire de Probabilites et Modeles Aleatoires CNRS-UMR 7599, Univ. Paris VI-VII, 4 
Pl.Jussieu, Paris, FRANCE 

E-mail address: ctoninel@ccr.jussieu.fr 



